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Abstract 

We have employed the semidiscrete variational generalized Peierls-Nabarro 
model to study the dislocation core properties of aluminum. The general- 
ized stacking fault energy surfaces entering the model are calculated by using 
first-principles Density Functional Theory (DFT) with pseudopotentials and 
the embedded atom method (EAM). Various core properties, including the 
core width, splitting behavior, energetics and Peierls stress for different dislo- 
cations have been investigated. The correlation between the core energetics 
and dislocation character has been explored. Our results reveal a simple rela- 
tionship between the Peierls stress and the ratio between the core width and 
atomic spacing. The dependence of the core properties on the two methods for 
calculating the total energy (DFT vs. EAM) has been examined. The EAM 
can give gross trends for various dislocation properties but fails to predict the 
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finer core structures, which in turn can affect the Peierls stress significantly 
(about one order of magnitude). 

§1. INTRODUCTION 

Dislocations which are one dimensional topological defects, are central to the understand- 
ing of mechanical properties of crystalline solids. The creation and motion of dislocations 
mediate the plastic response of a crystal to external stress. While continuum elasticity 
theory describes well the long-range elastic strain of a dislocation for length scales beyond 
a few lattice spacings, it breaks down near the singularity in the region surrounding the 
dislocation center, known as the dislocation core. The discrete nature of the real crystalline 
lattice avoids the conceptual difficulty posed by the continuum singularity and recovers 
the structural differentiation smoothed out by the continuum elasticity. There has been a 
great deal of interest in describing accurately the dislocation core structure on an atomic 
scale because of its important role in many phenomena of crystal plasticity (Duesbery and 
Richardson 1991, Vitek 1992). The core properties control, for instance, the mobility of 
dislocations, which accounts for the intrinsic ductility or brittleness of solids. The core is 
also important for the interaction of dislocations at close distances, which are relevant to 
plastic deformation. For example, by integrating the local rules derived from atomistic sim- 
ulations of core interactions into dislocation-dynamics simulations, a connection between 
micro-to-meso scales can be established to study dislocation reactions and crystal plasticity 
(Bulatov, Abraham, Kubin, Devincre and Yip 1998). 

Two types of theoretical approaches have been employed to study the core properties 
of dislocations. The first type is based on direct atomistic simulations employing either 
empirical potentials or first-principles calculations. Empirical interatomic potentials involve 
the fitting of parameters to a predetermined database and hence may not be reliable in 
describing the core properties, where severe distortions like bond breaking, bond formation 
and switching necessitate a quantum mechanical description of the electronic degrees of 
freedom. On the other hand, first-principles electronic structure calculations, though con- 
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siderably more accurate, are computationally expensive for studies of dislocation properties. 
The second type of approach is based on the framework of the Peierls-Nabarro (P-N) model 
which seems to be a plausible alternative to direct atomistic simulations. In fact, there 
has been a resurgence of interest in the simple and tractable P-N model for the study of 
dislocation core structure and mobility (Schoeck 1994, Joos, Ren and Duesbery 1994, Juan 
and Kaxiras, 1996). 

Peierls first proposed the remarkable hybrid model (1940) in which some of the details 
of the discrete dislocation core are incorporated into an essentially continuum framework. 
Nabarro (1947) and Eshelby (1949) further developed Peierls' model and gave the first 
meaningful estimate of the lattice friction to dislocation motion. Later attempts to generalize 
the original treatment of Peierls and Nabarro assumed a more general core configuration from 
which they derived the interactions between the glide planes which satisfy the Peierls integral 
equation. The essence of these models was captured in a more comprehensive approach by 
Vitek (1968, 1974), who introduced the concept of the generalized stacking fault: Consider 
a perfect crystal cut across a single plane into two parts which are then subjected to a 
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relative displacement through an arbitrary vector / and rejoined. The reconnected lattice 
will have a surplus energy per unit area 'y(f). As the vector f is varied to span a unit cell 
of the interface, ^y(f) generates the generalized stacking fault (GSF) energy surface. The 
procedure can be repeated for various crystal planes. The significance of the GSF surface 

— # 

(or 7-surface) is that for a fault vector / there is an interfacial restoring stress 

W) = -V( 7 (/5) (1) 

which has the same formal interpretation as the restoring stress in the P-N model. The 
P-N model has now come to represent a combination of the original continuum model and 
the GSF interplanar potential, and its accuracy can be affected by either component. At 
present, the GSF energies can be calculated using empirical interatomic potentials (like the 
embedded atom method - EAM) or electronic structure methods. While extremely useful 
as a conceptual framework, the P-N model becomes increasingly inaccurate for dislocations 
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with narrow cores, which is typically the case in covalently bonded solids (Joos, Ren and 
Duesbery 1994, Miller and Phillips 1996). The origin of this inaccuracy remains controver- 
sial, and it has not been unequivocally established whether or not the Peierls framework can 
be extended to capture such situations. Exploring the limits and extending the range of ap- 
plicability of the classic P-N models remains a worthwhile endeavor. Recently, a semidiscrete 
variational generalized P-N model has been proposed by Bulatov and Kaxiras (1997), which 
has been successfully implemented to a study of dislocation mobility in silicon. The new 
model predicts that the Peierls stress for the screw dislocation in Si is 0.065 eV/A 3 (using an 
energy surface from a classical interatomic potential), more than two orders of magnitude 
lower than the value 9.0 eV/A 3 obtained from the classic P-N model. A direct atomistic cal- 
culation using the same interatomic potential gives 0.021 eV/A 3 . The new model also gives 
satisfactory results for other core properties when compared to direct atomistic simulations 
for various dislocations in Si. 

The purpose of this paper is to apply this new model to aluminum which is a prototypical 
ductile metal with much lower Peierls energy and stress than silicon. On the other hand, 
just like silicon, aluminum is known to have a narrow core due to its large stacking fault 
energy. The successful application of the new model to aluminum will further prove its 
validity and versatility in predicting dislocation core properties of different materials. We 
have calculated the Peierls stress using both the semidiscrete generalized P-N model and 
atomistic simulations based on EAM. We have carried out systematic calculations of the 
core properties and the mobility of relevant dislocations in Al, and we have examined the 
relationship between the core properties (energetics, core width and Peierls stress) and the 
dislocation character, namely the angle between the dislocation line and its Burgers vector. 
In order to explore the dependence of the dislocation properties on the method employed 
for the GSF energy calculations, we have performed calculations using both first-principles 
electronic structure methods and the empirical EAM potential. 

The remainder of this paper is organized as follows: §2 describes the computational 
techniques used in our first-principles calculations for the GSF energy surface. §3 contains 
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a brief review of the new model on which this study is based. In §4 we present the results 
of the GSF energy surfaces for the (111) plane of aluminum using both the first-principles 
and empirical potential calculations. In §5 we compare the dislocation properties using the 
GSF surfaces from the two methods for obtaining the energetics. The correlation between 
the dislocation properties and the dislocation character is presented in §6, along with some 
general conclusions on the applicability of our approach. 

§2. COMPUTATIONAL METHODS 

The GSF energy surface is calculated within the framework of density functional theory 
(Hohenberg and Kohn 1964) in the local density density approximation (Kohn and Sham 
1965) to the exchange-correlation functional, using the expression proposed by Perdew and 
Zunger (1984). We refer to these calculations by DFT in the following. A kinetic energy 
cutoff of 12 Ry for the plane wave basis is used in the present calculations and the atomic 
structures are considered fully relaxed when the Hellmann-Feynman forces on each atom are 
smaller than 0.001 Ry/au. The calculated equilibrium lattice constant and bulk modulus 
are 3.94 A and 82.52 GPa, in good agreement with the corresponding experimental room- 
temperature lattice constant and bulk modulus of 4.05 A and 76.93 GPa, respectively. 

In addition to the lattice constant and bulk modulus, we have calculated the value of the 
intrinsic stacking fault energy and the unstable stacking fault energy. To simulate the block 
shearing process we employed a supercell consisting of six or nine atomic layers in the (111) 
direction. The intrinsic stacking fault configuration corresponds to a slip of a /\^6 in the 
(112) direction resulting in the stacking ABC|BCABC. The unstable stacking fault energy 
corresponds to the lowest energy barrier that needs to be crossed for the slip from the ideal 
configuration to the intrinsic stacking fault in the (112) direction. The calculations are per- 
formed at the theoretically determined in-plane lattice constant to eliminate any artificially 
imposed stress. For the reciprocal space integration we have used a k-point grid consisting of 
(16, 16, 4) divisions along the reciprocal lattice directions according to the Monkhorst-Pack 
scheme (Monkhorst and Pack 1976). This corresponds to 514 k-points in the entire Bril- 
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louin zone. Convergence tests were performed both for the number of divisions along each 
reciprocal space direction as well as the number of plane waves. Furthermore, the results for 
the 6-layer and 9-layer supercells are reasonably close, indicating adequate convergence with 
respect to the supercell size. Both atomic relaxations and volume relaxations were carried 
out to obtain accurate GSF energies. The value of the intrinsic stacking fault energy we 
obtained is 0.164 J/m 2 , in excellent agreement with the result of 0.165 J/m 2 of Sun and 
Kaxiras (1997), and that of 0.161 J/m 2 of Wright, Daw and Fong (1992). Experimental 
measurements range from a low value of 0.110 J/m 2 to a high value of 0.280 J/m 2 . Finally, 
the value of the relaxed unstable stacking fault energy obtained from our calculations is 
0.224 J/m 2 , which is identical to the value obtained by Sun and Kaxiras (1997). 

§3. SEMIDISCRETE VARIATIONAL PEIERLS-NABARRO MODEL 

To facilitate the presentation we adopt the following conventions: In Fig. 1, the XZ plane is 
the (111) glide plane, the Z axis is in the direction of the dislocation line, and the X axis is 
in the glide direction. The Y axis in normal to the glide plane. For planar dislocations, the 
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displacements along the Y direction are small. The Burgers vector b lies in the glide plane 
making an angle 8 with the Z axis. The Burgers vector is along the X axis (9 = 90°) for an 
edge dislocation and along the Z axis {6 = 0°) for a screw dislocation. The Burgers vector of 
a mixed dislocation has both an edge component, bsin9, and a screw component, bcos9. In 
general, the atomic displacements have components in all three directions rather than only 
along the direction of the Burgers vector. This is because the path along the Burgers vector 
may have to surmount a higher interplanar energy barrier in the GSF surface, that is, the 
GSF energy is reduced when the dislocation acquires additional displacement components 
in other directions. 

In the classic P-N formalism, the dislocation misfit is assumed to be confined into a 
single plane, the glide plane, separating two semi-infinite linear elastic continua. Between 
these two elastic half-spaces is placed a dislocation, conveniently represented as a continuous 
distribution of infinitesimal dislocations (Eshelby 1949) with density p(x). Here X is the co- 
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ordinate of the atomic row, which is always parallel to the dislocation line. A discrete lattice 
of arbitrary structure, deformed by the dislocation's displacement field, is superimposed on 
the elastic half-crystals. At a given point along the interface the resultant misfit, due to all 
the infinitesimal dislocations, is then balanced against the lattice restoring stress across the 
glide plane, F b (f(x)). This results in the P-N integrodifferential equation 
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Here, f(x) is the disregistry vector of the atomic row at point x related to the dislocation 
density by p(x) = df(x)/dx, and K is a constant depending on the elastic properties and 
the dislocation character. Its value for an isotropic solid is given by: 

u ,sin 2 9 9 . . 

where p and v are the shear modulus and Poisson's ratio, respectively, and 6 is the angle 
between the dislocation line and the Burgers vector. Thus, the pre-logarithmic elastic energy 
factor for a screw dislocation is K s = p/2ir, while for an edge dislocation it becomes K e = 
2n(i-v) • ^he dislocation density p(x) satisfies the normalization condition 

+00 +00 rif( >\ 

J p{x')dx' = J ^rdx' = b. (4) 

—00 —00 

If a simple sinusoidal form is assumed for F b (f(x)), as in the original P-N model, the 
disregistry vector is then given by the well-known analytical solution, 

/(*) = £ tan" 1 ^ + |, (5) 
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where 
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is the half width of the dislocation and F max is the maximum restoring stress. 

As pointed out recently by Bulatov and Kaxiras (1997), the classic P-N continuum model 
has the following flaws: (l)While the elastic energy between the infinitesimal dislocations is 



evaluated from a continuous integration, the misfit energy is sampled discretely across the 
glide plane in order to incorporate the discrete nature of lattice. Thus, the treatments of the 
two energy contributions are not on equal footing and the total energy is not variational. (2) 
The classic P-N model neglects the important degrees of freedom which participate actively 
in the translation of a dislocation over the Peierls barrier. (3) The elastic strain energy of a 
dislocation calculated within the P-N model can be unrealistically high, especially for solids 
with a narrow core. 

In order to resolve these problems, the semidiscrete variational generalized Peierls model 
has been developed by Bulatov and Kaxiras (1997). Within this approach the equilibrium 
structure of a dislocation is obtained by minimizing the dislocation energy functional 

U«* = E \^[Ke{^ P f + pf V? } ) + K sP f P f] 

+ £A* 73 (/1)-£^^(pN ) (7) 

i i,l 

with respect to the dislocation density or disregistry vector. Here, pf\ pf 1 and pf^ are the 
edge, vertical and screw components of the general interplanar displacement density at the 
i-th nodal point, and ^(fi) is the three-dimensional misfit potential. The corresponding 
applied stress components interacting with the pf\ pf^ and pf \ are = <t 2 i, = a 2 2 
and = (T23, respectively. K e and K s are the edge and screw pre- logarithmic energy 
factors defined earlier. The dislocation density at the i-th. nodal point is pi = (/j — — 
Xi-i), where fi and Xi are the disregistry vector and the coordinate of the i-th nodal point 
(atomic row), respectively. The remaining quantities entering in this expression are: Xij = 
\4>i,i-i(t>j,i-i + i'i-ij-i + ~ ~ with <j> itj =Xi- Xj, and ^ itj = In 

The first term in the energy functional, Eq.(7), represents the elastic energy, which has 
been discretized. Since any details of the displacements across the slip plane other than those 
on the atomic rows are disregarded, the dislocation density is constant between the nodal 
points. This explicit discretization of the elastic energy term removes the inconsistency in 
the original P-N model and allows the total energy functional to be variational. Another 
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modification in this approach is that the nonlinear misfit potential in the energy functional 
is a function of all three components of the nodal displacements, f{xj). Namely, in addition 
to the displacements along the Burgers vector, lateral and even vertical displacements across 
the slip plane are also included. This in turn allows the treatment of straight dislocations of 
arbitrary orientation in arbitrary glide planes. Furthermore, because the disregistry vector 
f(xi) is allowed to change during the process of dislocation translation, the Peierls energy 
barrier can be significantly lowered compared to the corresponding value from a rigid trans- 
lation. The response of the dislocation to an applied stress is achieved by minimization of the 
energy functional with respect to pi at a given value of the applied stress, rf\ An instability 
is reached when an optimal solution for p i no longer exists, which is manifested numerically 
by the failure of the minimization procedure to convergence. Within this formulation, we 
take as the definition of the Peierls stress the critical value of the external applied stress 
which produces this instability. 

§4. GENERALIZED STACKING FAULT ENERGIES 

The first-principles GSF energy surface 'j(f) for the (111) plane was calculated on a 
dense grid of 40 points in the irreducible part of the (111) slip plane (which is 1/12 of the 
area shown in Fig. 2). We have used a symmetrized polynomial basis to fit the GSF energy 
surface in order to evaluate the restoring force and facilitate the computation of dislocation 
properties. Since it is relatively much faster to calculate the GSF energy using the EAM 
method, we directly compute the GSF energy for any given disregistry vector. Although 
we do not need the EAM GSF energy surface explicitly for calculating any dislocation 
properties, we also fit that energy surface with the same basis of symmetrized polynomials 
for the purpose of comparison to the DFT calculation. The two energy surfaces are shown 
in Fig. 2. 

Both GSF surfaces maintain the rotation symmetry of the fee lattice. The three high 
peaks of the GSF surfaces correspond to the run-on stacking fault configuration ABC|CABC, 
in which the two C layers are neighboring to each other. The projection of the GSF energy 
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surfaces onto the [121] and [101] directions are shown in Fig. 3(a) and 3(b), respectively. 
Along the [121] direction, the first energy maximum encountered is the unstable stacking 
fault energy, which represents the lowest energy barrier for dislocation nucleation (Rice 
1992, Rice and Beltz 1994) and the first energy minimum at a /VQ (where a is the lattice 
constant) corresponds to the intrinsic stacking fault configuration, where a full dislocation 
dissociates into a pair of Shockley partials. On the other hand, the GSF surface projection 
along the [101] direction is symmetric with respect to the slip a /2^/2. For both the DFT 
and EAM calculations, the unstable stacking fault energy along [101] is found to be larger 
than that along [121]. This anisotropy of the unstable stacking fault energy will affect the 
emission of dislocations from a crack tip. The energy values for the various stacking faults 
obtained from the DFT and EAM calculations are summarized in Table I. The restoring 
stresses from both calculations and for different directions are shown in Fig. 4(a) and 4(b), 
respectively. 

§5. DISLOCATION PROPERTIES 

In this section we present the results of various dislocation properties for the screw, 30°, 

— # 

60° and edge dislocations with the Burgers vector, b = a/2 [101], using the GSF energy 
surface obtained from the DFT calculations and the direct EAM energy calculations. 

§5.1. Disregistry vector and dislocation width 

The disregistry vector / = fix + f 2 y + f3Z for the four dislocations can be calculated from 
the model using both the DFT and EAM energy surfaces (see Fig. 1 for the definition of X,Y 
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and Z). Since all four dislocations have the same Burgers vector b but different orientations 
(different angle of the dislocation line with respect to b), in order to examine the trend as a 

— * 

function of the angle one needs to project the disregistry vector / in two directions, parallel 

— * — * 

and perpendicular to the Burgers vector b. The components of / parallel and perpendicular 
to the Burgers vector are f a = f\ sin 9 + / 3 cos 9 and f b = —f\ cos 9 + / 3 sin 9, respectively. 
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The results for f a from both energy surfaces for the screw and 60° dislocations are presented 

— # 

in Fig. 5(a) and for the 30° and edge dislocations are shown in Fig 5(b) (all in units of |6|). 
In order to examine the degree of deviation of the disregistry vector f(f a ,f b ) from 

— * 

the direction of the Burgers vector b, we present in Fig. 6 the disregistry path for the four 
dislocations using the DFT (Fig. 6(a)) and EAM results (Fig. 6(b)). Although both methods 
predict that for the four dislocations the disregistry paths ) deviate significantly from 

— * — * 

the direction of b (vector AB in the figures) in order to lower the total energy, the magnitude 
of deviation obtained from the DFT calculation is much smaller than that from the EAM 
energy surface. More importantly, the deviation from the Burgers vector is almost zero in 
Fig. 6(a) for the atomic rows at the center of the dislocations whereas it reaches a maximum 
in Fig. 6(b). Apparently the difference between the two GSF energy surfaces results in the 
different displacement paths for the four dislocations. 

One can determine from Fig. 5 the dislocation half width (, defined as the atomic distance 
(nodal spacing) over which f a changes from (1/4)|6| to (3/4) The half width can also be 
calculated from the classical P-N model which assumes a sinusoidal form for the restoring 
stress F b (x) (Eq. (6)). Comparison of these two definitions provides information on how 
important are the details of the restoring stress in determining (. We have used a sinusoidal 
restoring stress with F max = F us , where F us is defined as the first maximum restoring stress 
encountered along the [121] direction for both the DFT and EAM cases. Values of the 
dislocation width of the four dislocations evaluated from the above two definitions are listed 
in Table II. We find that both DFT and EAM calculations give a dislocation core width 
that increases monotonically with the dislocation angle. The EAM calculation gives a wider 
dislocation core compared to that obtained from the DFT energy surface. This is due to 
the fact that the EAM potential gives smaller restoring stresses and GSF energies for a 
wide range of shears along the [121] direction, and in particular for the F us and the intrinsic 
stacking fault energy ^y is (see Fig. 3(a) and 4(a)). The smaller F us gives rise to a larger core 
width, while the smaller results in the dissociation of the full dislocation into partials 
and hence a wider core. 

11 



It is important to note that overall there is a good agreement (in particular for the 
DFT calculations) for the core width as evaluated from the two definitions described above. 
While the first definition takes into account the full details of the entire GSF surface in 
the evaluation of the disregistry vectors, the second definition involves only the maximum 
restoring stress F us . The agreement suggests that the details of the GSF surface are not 
important as far as the core width is concerned. The difference between the two definitions 
is larger for the 60° and edge dislocations when the EAM potential is used. This is because 
the EAM energy surface gives rise to a splitting into partials for both dislocations, which 
could only be manifested within the first definition. 

§5.2. Splitting into the partials 

The dislocation density for the screw, 30°, 60°, and edge dislocations calculated from 
the DFT GSF energy surface (solid curves) are presented in Fig. 7, together with the 
corresponding results from the EAM calculation (dashed and dotted curves). One can see 
that in all cases the DFT energy surface gives no splitting of the complete dislocations into 
partials. While the narrow double-peak structure found for the 60° complete dislocation is 
suggestive of a splitting, it is rather due to the fact that the nodal points along the X direction 
are not evenly spaced, i.e. they are distributed alternately by b p and b p /2, where b p is the 
Burgers vector of the Shockley partial. This in turn gives rise to density fluctuations over 
the neighboring atomic rows and hence a double-peak structure. A double-peak structure 
may indicate a splitting if the peaks are separated by a larger distance or the nodal points 
are evenly spaced, as in the cases of the 30° and edge dislocations. 

As far as these results are concerned, for the 30° and the 60° dislocations we plot both the 
edge and screw components. In contrast to the DFT results, the EAM calculations predict 
that the full edge and 60° dislocations will split into partials. The multi-peak found for the 
screw- dislocation suggests it might also be unstable. The 30° dislocation exhibits a narrow 
double-peak structure. The overall splitting trend found in the EAM calculations is due to 
the fact that the EAM potential gives a smaller intrinsic stacking fault energy. Realizing that 
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the EAM GSF energy surface differs from the DFT surface not only in the intrinsic stacking 
fault (ISF) energy, but also the energy profile around it, we have investigated the effect of 
the shape of the GSF surface on the splitting behavior for the full edge dislocation using the 
DFT GSF surface. We have kept the ISF energy of the DFT surface to be the real value 
of 0.164 J/m 2 , but varied the energy profile around the ISF to make it as flat as possible 
to be alike the EAM GSF surface (Fig. 8). We then calculated the dislocation density 
for the edge dislocation using the modified GSF surface, and we find that there is a small 
dissociation into partials (double-peak in the density distribution) with the spacing of one 
Burgers vector between the partials (Fig. 9). Therefore not only the ISF energy itself, but 
also the energy profile of the GSF energy surface around it could determine a splitting into 
partials. It is interesting to note that the EAM calculations predict a separation between 
partials of about 8-9 A (about three times the Burgers vector), in agreement with the values 
of about 6-7 A predicted by the elastic continuum theory (Hirth and Lothe 1982). It is also 
interesting to examine the character of the resultant partials. The complete edge dislocation 
dissociates into two symmetric 60° partials, whereas the 60° dislocation dissociates into a 
30° (left double-peak) and a 90° (right double-peak) partial. The double-peak structure 
again stems from the inequivalent nodal spacing between neighboring atomic planes. The 
fact that the density of the screw component (dotted curve) vanishes at the point where the 
edge component (right double peak) reaches its maximum, indicates the pure edge character 
of the 90° partial. The 30° dislocation exhibits a weak tendency for dissociation into a screw 
and a 60° partial, while the edge component for the 60° partial is singly peaked. Our results 
for the dissociation of the full 60° dislocation into partials are in agreement with those from 
the direct atomistic simulations by Bulatov (1998) using the same EAM potential. 

§5.3. Energetics and lattice resistance 

In Table II we present the results for the the core energy, defined as the sum of the elastic 
and misfit energies, the separate contributions to the core energy from the elastic and GSF 
energies, and the Peierls stress for the four dislocations using the DFT and EAM potentials. 
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Also for comparison, we include in Table II the values of the Peierls stress for the screw 
and 60° dislocations obtained from the direct atomistic simulations using the same EAM 
potential. 

The results for the energetics from the DFT calculations are overall in good agreement 
with those from the EAM calculations, including the trend across the types of dislocations. 
We identify the configuration dependent part of the total energy as the core energy, including 
the elastic energy and GSF energy. The GSF energy contribution (positive in sign) increases 
monotonically in going from the screw to the edge dislocation; the elastic contribution (neg- 
ative in sign) decreases as the angle of the Burgers vector with respect to the dislocation line 
increases. The elastic energy, ignored in some previous studies, is the dominant contribution 
to the core energy (about a factor of two larger than the GSF energy), and more importantly 
it is strongly dependent on the dislocation character. 

The mobility of dislocations is characterized by the Peierls stress and the Peierls energy, 
the former being defined usually as the maximum derivative of the latter with respect to 
dislocation translation. As mentioned earlier, the Peierls stress in this work is obtained 
by evaluating the critical value of the applied stress r, at which the dislocation energy 
functional fails to be minimized with respect to pi through standard conjugate gradient 
techniques. This approach is more accurate and physically transparent because it captures 
the nature of the Peierls stress as the stress at which the displacement field of the dislocation 
undergoes a discontinuous transition. 

The results for the Peierls stress, listed in Table II, calculated from the DFT GSF energy 
surface are in good agreement with those from the direct EAM atomistic simulations and 
those from the EAM calculations except for the edge dislocation. For the edge dislocation, 
the EAM calculations yield a Peierls stress which is an order of magnitude larger than 
that from the DFT calculations. This discrepancy could be due to the dissociation of 
the full edge dislocation into partials. The effects of the dissociation on the Peierls stress 
have been studied by Benoit et al. (1987). These authors argued that if the equilibrium 
separation between the two Shockley partials is an integer or half integer multiple of the 
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vector (a/2) (110), so that the two partials move in phase and reach the troughs or the 
crests of the Peierls potential simultaneously, then the Peierls stress required to move the 
extended configuration is simply the stress required to move an isolated partial dislocation. 
On the other hand, if the two partials are rigidly separated by (1/4) or (3/4) times the 
vector (a/2) (110), then the partials are exactly out of phase and the Peierls stress on them 
are always equal and opposite, and hence the applied stress required to move the rigid 
configuration vanishes. 

In order to understand the discrepancy for the Peierls stress of the edge dislocation from 
the two calculations, we next investigate the effects of the splitting on the Peierls stress. 
Since the EAM gives a dissociation for the perfect edge dislocation into two symmetric 60° 
partials with a separation of about three times the Burgers vector (see §5.2), one can infer 
from the work of Benoit et al. that the Peierls stress for each of the two partials is equal to 
the Peierls stress required to move the extended dislocation, i.e. the two partials and the 
intrinsic stacking fault in between. Continuum elastic theory gives that the separation of 
the partials is 

d=^ 2 -±^, (8) 

where b p is the partial Burgers vector and 7 is is the intrinsic stacking fault energy. Since 
the DFT GSF energy surface predicts no splitting for the edge dislocation, one can force the 
dissociation into partials by either increasing the elastic constant K in Eq. (3) or reducing 
the intrinsic stacking fault energy. 

We find that increasing the pre-logarithmic factor K in the elastic energy term (by up 
to a factor of 10) but keeping the same DFT GSF surface gives rise to a wider dislocation 
density distribution but yields no dissociation. On the other hand, if we keep the same 
pre-logarithmic factor K but reduce the DFT value for the intrinsic stacking fault 7 is from 
0.164J/m 2 to 0.096 J/m 2 (the vicinity of the stacking fault in the GSF surface is also 
reduced by this rescaling of the energy), we find that the perfect edge dislocation splits into 
two 60° partials separated by one Burgers vector and a Peierls stress of 0.33 meV/A 3 , an 

15 



order of magnitude larger than the original value of 0.02 meV/A 3 . If we further reduce 
the stacking fault energy to 0.085 J/m 2 , we find that the edge dislocation splits into two 
partials separated by a distance of 36, but with the same Peierls stress (0.33 meV/A 3 ) for 
the extended configuration. Since our calculations give separations between partials that 
are integer multiples of the Burgers vector, we can infer from the work of Benoit et al. that 
the stress of 0.33 meV/A 3 should be equal to the Peierls stress for moving an isolated 60° 
partial. Thus, we find that the origin of the discrepancy of the Peierls stress for the edge 
dislocation between the DFT and EAM calculations is that the DFT calculations predict 
no splitting for the edge dislocation and hence a very low Peierls stress, whereas the EAM 
calculations predict a splitting into partials due to the smaller stacking fault energy, and 
hence a much higher Peierls stress. It was not possible to test whether the Peierls stress 
for an isolated 60° partial is equal to the Peierls stress for the extended dislocation (0.33 
meV/A 3 ), because an isolated partial dislocation introduces an infinite stacking fault on the 
glide plane. 

In order to test the accuracy of the present model for predicting the Peierls stress, we 
also have performed direct atomistic simulations using the same EAM potential for both 
the screw and 60° dislocations. The values for the Peierls stress from the direct atomistic 
calculations listed in Table II are in very good agreement with our model calculations. For 
comparison, we have also calculated the Peierls stress for the screw and 60° dislocations 
using the expression of Joos and Duesbery (1997) 

a p = -r ex py — ir)> ( 9 ) 

(Jj (Jj 

where a' is the atomic spacing, while the core width ( is determined from the EAM calcu- 
lations. Eq. (9) gives a Peierls stress 1.44xl0 -2 meV/A 3 for the screw dislocation, and of 
4.8xl0~ 4 meV/A 3 for the 60° dislocation, both several orders of magnitude smaller than the 
corresponding values found from the current model and the results of the direct atomistic 
simulations. 

§6. CORRELATION BETWEEN DISLOCATION PROPERTIES AND 
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DISLOCATION CHARACTER 



Using the DFT GSF surface we have studied the dislocation properties of nineteen different 
dislocations that have the same Burgers vector but different orientations, in order to examine 
the possible correlation between the dislocation properties and the dislocation character. 
The angle between the dislocation line and the Burgers vector has been varied from 0° to 
90°. 

In Fig. 10 we present the dislocation core width as a function of the dislocation angle 
using the two definitions described in §5.1. The width increases monotonically with the 
dislocation angle. Even though the first definition takes into account the details of the 
entire GSF surface while the second definition employs only a single special point of a 
simpler sinusoidal GSF (namely F us ), the overall agreement is very good (less than 10% 
difference) . 

The core energy, along with its separate energy contributions from the elastic energy 
(interaction between the infinitesimal parallel dislocations on the glide plane) and the GSF 
energy are presented in Fig. 11 as a function of the angle. We find that the core energy 
and the elastic energy decrease monotonically as the angle increases from 0° (screw) to 
90° (edge), whereas the GSF energy increases with the angle. As the angle increases the 
pre-logarithmic factor K of a mixed dislocation increases and hence the elastic energy is 
lowered. On the other hand, the monotonic increase of the GSF energy with respect to 
the angle is due to fact that the core width increases monotonically with respect to the 
dislocation character. The elastic energy dominates (in absolute value) the GSF energy 
and increases in magnitude faster with increasing angle than the GSF energy, indicating a 
weaker anisotropy for the GSF energy compared to that for the elastic energy. Thus, the 
elastic energy not only is the dominant contribution to the total energy stored in the core 
region, but also is more sensitive to the dislocation character than the GSF energy. 

Finally, we have calculated the Peierls stresses for all the nineteen dislocations in an 
effort to correlate the Peierls stress with the dislocation character. Most of the dislocations 
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in the fee lattice have non-even nodal spacings, except for the 30° and edge dislocations. In 
Fig. 12 we plot \n(o p a/ Kb) as a function of (/a, where ( is the half width of the dislocation 
core calculated from the atomic spacing over which the disregistry vector changes from \b\/4 
to (3/4) | b\, and a is the average nodal spacing along the X direction. We also show for 
comparison the analytic expression Eq. (9) of Joos and Duesbery (1997), with a dashed line 
(slope of -2ir). It should be emphasized that this expression is valid only for evenly spaced 
atomic planes. Most of the calculated values can be fitted with 

<y v = —e-^i\ (10) 

(X 

shown as the solid line. The points corresponding to the 30° and edge dislocations deviate 
from the above line because of the even nodal spacing. We have recently formulated a 
general theory to account for the non-evenly spaced dislocations (Lu, Kioussis, Bulatov and 
Kaxiras, 1999) which shows that the Peierls stress for evenly-spaced dislocations is lowered 
by several orders of magnitude compared to that of the non-evenly spaced dislocations. 
On the other hand, the deviation from the common trend of the values for the 10.9° and 
14.9° dislocations is unclear at present. The Peierls stress is more sensitive to the average 
atomic spacing a than to the half width. For example, while both the screw and 14.9° 
dislocations have predominant screw components and similar half widths of 2.1 A and 2.3 
A , respectively, they have quite different atomic spacings, 1.2 A and 0.3 A , respectively. 
This results in a Peierls stress of 0.04 meV/A 3 for the 14.9° dislocation, almost two orders 
of magnitude smaller than that of 1.60 meV/A 3 for the screw dislocation. 

In conclusion, we have performed DFT and EAM calculations to obtain the GSF en- 
ergy surfaces for the (111) glide plane of Al. From those calculations we extracted the 
properties for various dislocations, using the semidiscrete variational generalization of the 
Peierls-Nabarro model. We have demonstrated that although the EAM gives the general 
trends for various dislocation properties, it fails to predict the finer structure of the dislo- 
cation core, such as the dissociation into partials, which in turn determines the mobility of 
dislocations. Since the splitting of a full dislocation into partials depends strongly on the 
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intrinsic stacking fault energy, direct atomistic simulations based on empirical potentials 
may also fail to predict the correct dissociation behavior. Accordingly, the results of the 
present work indicate that the accurate DFT calculations of the GSF surface combined with 
the semidiscrete variational P-N model enable the study of dislocation core properties more 
accurately and expediently. Moreover, this model could be extended to study a wide range 
of problems that are associated with more complex dislocations processes such as cross slip, 
dislocation intersection, reaction, etc. 
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TABLES 

TABLE I. The energetics and fault vectors for the four stacking faults obtained from the DFT 



and EAM calculations. All energies in J/m 2 . 

Vector DFT EAM 

Intrinsic stacking 1/6[121] 0.164 0.120 

Unstable stacking 1/10[121] 0.224 0.141 

Unstable stacking 1/4[101] 0.250 0.663 

Run-on stacking 1/3[121] 0.400 1.354 
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TABLE II. Core widths ( (in A ), with £ 7 from the definition through the disregistry vector 
and ( n from the definition through the restoring stress; total core energies and separate contri- 
butions from the elastic and GSF energies (in eV/A); and Peierls stress (in meV/A 3 ), for the four 
dislocations from both the DFT and EAM calculations. The last two rows list the Peierls stress 
for the screw and 60° dislocations from the direct atomistic simulations (At. Sim.) using the same 
EAM potential and from the expression of Joos and Duesbery (JD). 







screw 


30° 


60° 


edge 


Core widths (C 7 ,C 7/ ) 


DFT 


(2.1,2.1) 


(2.5,2.4) 


(3.0,2.9) 


(3.5,3.2) 




EAM 


(3.7,3.4) 


(4.3,3.8) 


(5.4,4.7) 


(6.4,5.2) 


Core energy 


DFT 


-0.0834 


-0.1096 


-0.1678 


-0.1979 




EAM 


-0.0534 


-0.0813 


-0.1413 


-0.1730 


Elastic energy 


DFT 


-0.1828 


-0.2317 


-0.3199 


-0.3666 




EAM 


-0.1679 


-0.2086 


-0.2960 


-0.3445 


GSF energy 


DFT 


0.0938 


0.1221 


0.1521 


0.1688 




EAM 


0.1145 


0.1273 


0.1547 


0.1716 


Peierls stress 


DFT 


1.60 


0.33 


0.61 


0.02 




EAM 


0.55 


0.21 


0.28 


0.15 




At. Sim. 


0.51 




0.29 






JD 


0.0144 




0.00048 
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FIGURES 



FIG. 1. A Cartesian set of coordinates showing the directions relevant for dislocations in Al. 

FIG. 2. The GSF energy surfaces for displacements along a (111) plane in Al (J/m 2 ) (the 
corners of the plane and its center correspond to identical equilibrium configurations, i.e., the ideal 
Al lattice): (a) from DFT pseudopotential plane- wave calculations; (b) from EAM calculations. 
The two energy surfaces are displayed in exactly the same perspective and on the same energy 
scale to facilitate comparison of important features. The EAM energy surface extends to values 
about 4 times larger than the maximum value displayed (black truncated regions). 

FIG. 3. The projections of the GSF energy surfaces onto (a) [121] and (b) [101] directions from 
both the DFT and EAM calculations. 

FIG. 4. The restoring stresses along (a) [121] and (b) [101] (b) directions calculated from both 
the DFT and EAM GSF energy surfaces. 

FIG. 5. The disregistry vectors in unit of the Burgers vector obtained from the DFT and EAM 
calculations for: (a) the screw and 60° dislocations; (b) the 30° and edge dislocations. 

FIG. 6. The disregistry paths for the four dislocations obtained from (a) DFT calculations; (b) 
EAM calculations. 

FIG. 7. The density distributions for the four dislocations (clockwise): screw, 30°, 60° and 
edge, obtained by using DFT and EAM calculations: For the 30°, 60° and edge dislocations, the 
characters of resultant partials within the EAM calculations are indicated (see also text). 

FIG. 8. The projection of the DFT GSF energy surface along [121] direction with and without 
modifying the real energy profile. The solid curve represents the real energy profile without modi- 
fication while the dashed curve represents the modified GSF energy surface but keeping the same 
intrinsic stacking fault energy as its real value. 
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FIG. 9. The dislocation density distribution for the edge dislocation calculated from the modi- 
fied GSF energy surface (dashed line) and the real GSF energy surface without modification (solid 
line). The distance is in terms of the full Burgers vector. 

FIG. 10. The half width (in A) as a function of dislocation orientation using the two different 
definitions. The square stands for the definition through the disregistry vector and the dot stands 
for the definition through the restoring stress (only DFT results shown). 

FIG. 11. The core energy, elastic energy and GSF energy as a function of dislocation orienta- 
tions. 

FIG. 12. The scaled Peierls stress as a function of the core width and the atomic spacing 
perpendicular to the dislocation line (the dashed line corresponds to the formula from Joos et al). 
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